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We discuss the influence of disorder on domain structure formation in ferromagnetic phase of di- 
luted magnetic semiconductors (DMS) of p-type. Using analytical arguments we show the existence 
of critical ratio v cr of concentration of charge carriers and magnetic ions such that sample critical 
thickness L CT (such that at L < L CI a sample is monodomain) diverges as v — > u CI . At v > v CT the 
sample is monodomain. This feature makes DMS different from conventional ordered magnets as 
it gives a possibility to control the sample critical thickness and emerging domain structure period 
by variation of v. As concentration of magnetic impurities grows, v CT — > oo restoring conventional 
behavior of ordered magnets. Above facts have been revealed by examination of the temperature of 
transition to inhomogeneous magnetic state (stripe domain structure) in a slab of finite thickness L 
of p-type DMS. Our analysis is carried out on the base of homogeneous exchange part of magnetic 
free energy of DMS calculated by us earlier [Phys. Rev. B , 67, 195203 (2003)]. 

PACS numbers: 72.20.Ht,85.60.Dw,42.65.Pc,78.66.-w 



The structure of domains and domain walls in conven- 
tional ordered magnets has been well studied both ex- 
perimentally and theoretically several decades ago (see, 
[l|,|2| and references therein). On the other hand, the im- 
portant question about the domain structure formation 
in the ferromagnetic phase of p-doped diluted magnetic 
semiconductors (DMS) 0, (which can be regarded as 
disordered magnets) has not been addressed yet. The 
characteristics of the domain structure in the DMS films 
of III-V type have been investigated in Ref (see also 
references therein). In our opinion, the main problem 
was the lack of suitable "continuous" free energy func- 
tional of DMS, which is necessary to theoretically study 
their domain structure. Recently [jj such free energy 
functions have been derived microscopically from Ising 
and Heisenberg models for DMS. In the language of phc- 
nomenological theory of magnetism these functions corre- 
spond to so-called homogeneous exchange parts of total 
phenomenological free energy of DMS. To describe the 
domain structure properties, these contributions should 
be completed by inhomogeneous exchange and magnetic 
anisotropy energies. It was demonstrated experimentally 
(see |(| and references therein) that magnetic anisotropy 
exists in DMS of (Ga,Mn)As type. At the same time it 
was demonstrated in jfj that unstrained samples (which 
can be well described by Heisenberg model Q) have easy 
plane magnetic anisotropy, while uniaxially strained sam- 
ples (Ising model 0) have anisotropy of easy axis type. 

It is well known (see, e.g. that at low tempera- 
tures the domain pattern formation is primarily due to 
the rotation of magnetization vector with constant modu- 
lus being saturation magnetization Mq. On the contrary, 
for the temperatures close to T c , this structure is formed 
by the variation of modulus of M rather then its rota- 
tion. This means that above homogeneous exchange part 
of the magnetic energy of DMS will only contribute to its 



domain structure in the vicinity of T c . At low temper- 
atures the influence of disorder on domain structure of 
DMS will be small so that it will resemble very much the 
domain structure of conventional magnetically ordered 
substances. 

In the present paper we suggest a theory of inhomo- 
geneous magnetic state (stripe domain structure) in the 
DMS slab in the vicinity of ferromagnetic phase transi- 
tion temperature. We analyze the sample of finite thick- 
ness L. We show that the impurity character of ferro- 
magnetism in DMS results in substantial narrowing of 
the region of temperatures and sample thicknesses, when 
domain structure exists. For example, beyond mean field 
approximation (i.e. when disorder in magnetic ions sub- 
system becomes substantial) even at zero temperature 
domain structure appears not at L = (as in the case of 
ordered magnet, see 0,0), but at some threshold value 
L iT , depending on the ratio v of charge carriers and mag- 
netic ions concentrations, see Fig.l. This effect makes the 
domain structure of DMS (disordered magnets) qualita- 
tively different from that of conventional ordered mag- 
nets. The developed formalism can be easily applied for 
thin DMS films. 

Consider the slab-shaped sample of DMS with slab 
thickness L (Fig.l, lower panel). Let z axis is magnetic 
anisotropy axis (and xy plane is the plane of anisotropy 
for Heisenberg model). The phenomenological free en- 
ergy of DMS near T c can be written in the form (see, e.g. 

0) 

F = Jdvl~a (VM) 2 + /£ N + f p (M) - ±MH D } , 

where M is a magnetization vector, a is inhomogeneous 
exchange constant, H jj is a demagnetizing field, P stands 
for H (Heisenberg model) or / (Ising model), f p N are 
the anisotropy energies and f p (M) are the homogeneous 
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exchange energies. For Heisenberg model (easy plane 
anisotropy) 
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where m = \m\, m = < S >/S, S* is a spin of a magnetic 
ion. Bar means the averaging over spatial disorder in 
magnetic ion subsystem in DMS .while angular brackets 
mean the thermal averaging, see ||| for details. The rela- 
tion between m and M in this case is usual: m=M/Mo, 

where Mq = < S >(T = 0) is saturation magnetization. 
For Ising model (easy axis anisotropy) 



/an = \(i {M 2 X + M 2 y ) , (3 > 0, 
f I (m) = ±m i (l-Ai) + ±m*A{ 
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^0,1 are the real and imaginary parts of Fourier image of 
distribution function of random magnetic fields, acting 
among magnetic impurities in DMS, J(r) is an interac- 
tion between spins of magnetic impurities in the Heisen- 
berg or Ising Hamiltonians (this is actually RKKY inter- 
action, see 00 for details), n(r) is (spatially inhomoge- 
neous) concentration of magnetic ions. Index 1/2 in the 
functions © means that free energies and |JB} have 
been derived for spin 1/2 (see [7j). However, our analysis 
shows that all our results remain qualitatively the same 
for arbitrary spin. 

Equilibrium distribution of magnetization in DMS can 
be obtained from the equation of state SF/Sm — and 

Maxwell equations rot h — 0, div (ji, + Anrfij = with 
boundary conditions for slab geometry 
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where h = Hp, /Mo, is a demagnetizing field in a vac- 
uum. It was shown in Q that for sufficiently thick slabs 



and P < 4-7T the following equation for distribution 
of magnetization in DMS can be derived from the above 
equation of state with respect to JD 
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where [x± = 1 + 4tt//3, b H = 1 - 2Af , c H = Ag /5, 
bj = 1— A\, ci = ylg / 6, where in the ferromagnetic phase 
(T < T c ) functions b P (T) < (P = H, I). It should be 
noted here, that the different forms of anisotropy ener- 
gies for Heisenberg and Ising (JSJ models in the above 
suppositions do not influence the form of equation 1)110. 

It was shown in (see also 0), that transition from 
paramagnetic phase to the ferromagnetic phase with do- 
main structure (domain state) occurs via phase transition 
of second kind. This means that for the determination of 
transition temperature Tk to the domain state it is suffi- 
cient to consider the linearized version of Eq. (|llfl . We 
look for its solution in the form 



m z = Acosqzcoskx. 



(12) 



Here cos qz determines the spatial inhomogeneity along 
z direction, while coskx defines a "linear" domain struc- 
ture (in x direction so that domain walls lie in yz plane) 
with a period d = 2ir/k. 

Substitution of solution (|12)l into the linearized version 
of (|11|) gives the equation relating q and k 
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Cp = -h 
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To obtain the dependence of Tk on sample thickness, we 
need one more equation relating £ (and by its virtue 
Tk, see (01,(0), k and q. Such equation can be got- 
ten substituting 1)12(1 into boundary conditions 110(1 with 
respect to vacuum solutions hi^ = C exp (— k\z\) cos kx, 
hx = C exp (— k\z\) sin kx. It reads 
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Equations ((13(1 and 1141) constitute a close set of equa- 
tions for the instability temperature Tk and equilibrium 
domain structure period as a function of sample thickness 
L and concentration ratio v. The equations 113() and 1(14(1 
can be reduced to a single equation for the dependence 
of bp on y — ky/a 
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(Cp - V 2 ) = ™ + - arctan ^ — — , 
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where rj 2 = L 2 /a. The equation ((15(1 is a single equation 
for C(y) at different r\. This dependence has the form of 
a curve with a minimum. The real transition to domain 
state in DMS occurs when Tk reaches its maximal value 
as a function of y. This can be demonstrated by sub- 
stitution of a solution of nonlinear equation 1(11(1 in the 
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form of infinite series in small parameter proportional to 
\T— Tfc| into free energy Q with its subsequent minimiza- 
tion over y similarly to Ref.0. Since (p = — bp for both 
models is a decreasing function of temperature (e.g. for 
both models in a mean field approximation ( = 1/t — 1 , 
r = T/T cMF 0), the coordinates of minimum Cp m 
and y min of the curve determine the equilibrium 
temperature of a phase transition to the domain state 
X;- and equilibrium period of emerging domain structure 
A = 2ir/y min as functions of dimensionless sample thick- 
ness r\. 
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FIG. 1: Upper panel - schematic phase diagram of the DMS 
slab in the coordinates temperature - thickness. Region 1 cor- 
responds to paramagnetic phase for both mean field (MF) and 
random field (RF) approximations, region 2 - paramagnetic 
phase for MF and domain state for RF, region 3 - domain 
state for both cases. The nonzero threshold sample thick- 
ness occurs for RF case. Horizontal asymptotes correspond 
to T {y) — T C /T C MF, T c is a temperature of phase transi- 
tion into ferromagnetic homogeneous (i.e. domainless) state. 
Lower panel - geometry of the sample. 

We have taken the minimum of implicit function 
C(v) C3l numerically to get the dependencies A(ry) and 
Cp la (v)- They are reported on inset to Fig. 2. ft is seen 
that dependence Cp ln (rj) decays rapidly as r\ — ► oo, at 
T) — > Cp lm — * oo- It can be shown that at large ij 
A oc 



y/rj = sjLl^fa. This behavior is typical for stripe 
domain structure in ordered magnets (see @>E|) an d is 
seen on the figure. 

Now we consider the explicit dependencies bp{y,r), 
v = n c /m = x c /xi, x c j — n Cj if2, Q is DMS unit cell 
volume, n c and rii are charge carriers (electrons or holes) 
and magnetic ions concentrations respectively. We have 
in dimensionless variables 



*(f ) = cxp 
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p = h,i, 



where £ = J xt /3 t/2T. In the expressions (|16|l we used 
RKKY potential in the simplest possible form corre- 
sponding to one-band carrier structure 



J(f) = -J x* /3 R(2k f r), F(x) 



x cos x — sin x 



, (17) 



where J Q = (3/tt) 1/3 (3/27i 2 ) J^n 2 / 3 m d , J ci is a carrier- 
ion exchange constant, rrid is the density of states effec- 
tive mass. Functions ¥>o(£) and </?i(£) have the form 



<A) (0 + ifi (0 



/ {l-exp(-t^F(y))}y 2 dy. 
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The dependencies Ch{t) and Cr( r ) a t different v are 
shown on Fig. 2. In a mean field (MF) approximation 
(hmf = (imf — Cmf = 1/t — 1 is unbounded at T = 0, 
while beyond this approximation functions Qh and Q 
have finite values at T = 
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FIG. 2: Left panel: dependence of £h and £j on normalized 
temperature at different (A Points where £p = correspond to 
T c (f). Dashed line labelled MF corresponds to £mf = 1/t — 1. 
Right panel: equilibrium period of emerging domain struc- 
ture and parameter ( versus dimensionless sample thickness 
at n± — 10 . 

These finite values are indeed "responsible" for the 
emergence of threshold sample thickness L tr beyond 
mean field approximation. Having dependencies l|16(l . we 
can solve them numerically for Cp m (determined above 
from Q15[0 to obtain the dimensionless phase transition 
temperature = Tk /T c m f as a function of critical sam- 
ple thickness L cr -/^/a. In other words, here we have the 
phase diagram of DMS slab in the coordinates (T,L). 
This is reported on Fig. 3. 

The presence of L CT (v) is clearly seen. The asymptotes 
for large L cr are due to the dependence of equilibrium 
(i.e. to the ferromagnetic phase without domain struc- 
ture) phase transition temperature r c = T c /T c mf on the 
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FIG. 3: Phase diagram of DMS slab in the coordinates T, L. 
Horizontal asymptotes - T c (v) similar to Fig.l, fi± — 10. 



concentration ratio v. Latter dependence is given by the 
conditions Qh — and = (see ©, ®) for Heisen- 
bcrg and Ising models respectively. It was shown in Ref. 
8] that impurity ferromagnetism in DMS is possible for 
< v < v? x (P = H, I, = 0.0989, i£ = 0.2473 @) so 
that t c {v ct ) = 0. This means that t c decays as v grows 
and the region (T,L), where ferromagnetic domain state 
exists in DMS, diminishes substantially (compared to the 
case of ordered ferromagnets) and vanishes as v — > v CT . 




FIG. 4: Threshold sample thickness (i.e. that at T — 0) and 
corresponding domain structure period as functions of con- 
centration ratio v at fi± = 10. Critical ratios for Heisenberg 
and Ising vl r are shown. For mean field approximation 
(ordered magnet) u C i — > oo. 



This fact makes it possible to control the critical thick- 
ness of DMS sample by changing v. This, in turn, might 
give a possibility to engineer the domain structure in 
nanocrystals od DMS, which is useful for many technical 
applications (see 0,0 and references therein). Note that 
our formalism permits to calculate v CT and other charac- 
teristics of domain structure for any temperature (not 
only neat T c and sample geometry). 

Now we present some numerical estimations. The ma- 
jor problem here is uncertain value of inhomogeneous 
exchange constant a. It can be estimated by the ex- 
pression a w kBT c a 2 / (M s fio), where a = AA is a 
typical value of lattice constant for DMS, M s « 50 mT 
|6| is a saturation magnetization (of localized spin mo- 
ments) of DMS, T c = T cM F ~ 100Jf 3 is a temperature 
of transition to homogeneous ferromagnetic state in a 
mean field approximation, fee and /io are Boltzmann con- 
stant and Bohr magneton respectively. Evaluation gives 
a 1 / 2 200A. From Fig. 4 for v — 0.075 we have thresh- 
old sample thickness L tT ~ 50a 1 / 2 = lOOOOA = \[xm and 
corresponding period of domain structure is ~ 25a 1 / 2 = 
5000A = 0.5^m for Heisenberg model. The same values 
for Ising model occur at v ~ 0.2. These values are in fair 
agreement with results of Ref. 0- Moreover, for differ- 
ent v we have quite different values of Ltr and y tT . This is 
the base for above discussed domain structure engineer- 
ing. For our picture to give the quantitative description 
of experiment in real DMS, the precise experimental de- 
termination of inhomogeneous exchange constant a and 
anisotropy constant /3 is highly desirable. 

Here we presented a formalism for the calculation of 
properties of domain structure in DMS. Our present re- 
sults about phase diagram of DMS is the simplest appli- 
cation of the formalism. Generally, it permits to calcu- 
late all desired properties of domain structure (like the 
temperature and concentration dependencies of domain 
structure period and domain walls thickness) in the entire 
temperature range as well as to account for more complex 
then slab sample geometries. Latter can be accomplished 
by applying different from (|10|l boundary conditions. The 
external magnetic field can also be easily taken into ac- 
count. However, far from the solution of resulting 
nonlinear differential equations would require numerical 
methods. 



Note, that in MF approximation it is very easy to solve 
(|TT))l analytically to get t c mf = V(Cp + !)■ Resolving 
the equations (|18fl for £p m we obtain the dependence of 
threshold thickness on concentration ratio v. This depen- 
dence along with corresponding period of domain struc- 
ture is shown on Fig. 4 for Heisenberg and Ising models. 
It is seen that as v approaches u cx Ltr — * oo so that DMS 
sample looses its ferromagnetism (both homogeneous and 
inhomogeneous) . The period of domain structure also di- 
verges as v — > i/ CT . 
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